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Abstract 

In the series of models with interacting particles in stochastic geometry, a further 
contribution presents the facet process which is defined in arbitrary Euclidean 
dimension. In 2D, 3D specially it is a process of interacting segments, flat sur¬ 
faces, respectively. Its investigation is based on the theory of functionals of finite 
spatial point processes given by a density with respect to a Poisson process. The 
methodology based on L 2 expansion of the covariance of functionals of Poisson 
process is developed for 17-statistics of facet intersections which are building 
blocks of the model. The importance of the concept of correlation functions of 
arbitrary order is emphasized. Some basic properties of facet processes, such 
as local stability and repulsivness are shown and a standard simulation algo¬ 
rithm mentioned. Further the situation when the intensity of the process tends 
to infinity is studied. In the case of Poisson processes a central limit theorem 
follows from recent results of Wiener-Ito theory. In the case of non-Poisson pro¬ 
cesses we restrict to models with finitely many orientations. Detailed analysis 
of correlation functions exhibits various asymptotics for different combination 
of 17-statistics and submodels of the facet process. 
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1 Introduction 

Modeling of particle systems with interactions by means of finite point processes 
presents an interesting field of stochastic geometry the research of which is not 
yet completed. The work of [6] concerning area-interaction models of planar 
discs was later developed by [9] to models of interacting discs with more gen¬ 
eral densities (with respect to the Poisson process) from exponential families. 
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The authors also developed sophisticated simulation procedures of the models. 
The motivation for the modeling was to serve for statistical purposes in the 
evaluation of real data. 

In the present paper we develop similar interaction processes in a bounded 
window in arbitrary Euclidean dimension. The difference in comparison to fiJj 
is that the particles (called facets) are lower dimensional, they form compact 
subsets of hyperplanes. The interactions arise from intersections of facets and 
the their global amounts form [/-statistics [T2] - In 2D and 3D these models 
may also serve to real data evaluation of segment, surface processes, respectively. 
Rather than statistics we develop the facet process theory in arbitrary dimension 
and deal with limit behavior when the intensity of the reference Poisson process 
tends to infinity. 

Recently in [2] functionals of spatial point processes given by a density with 
respect to the Poisson process were investigated using the L 2 expansion from [8] 
which is applied to the product of a functional and the density. Using a special 
class of functionals called [/-statistics closed formulas for mixed moments of 
functionals are obtained. Similar formulas, but under a stronger assumption 
of a product form of the driving function of the functional, were derived in [3] 
using the Georgii-Nguyen-Zessin formula. In processes with densities the key 
characteristics are the correlation functions |3] of arbitrary order which are dual 
to kernel functions of the density as a function of the Poisson process. 

We call facets some bounded subsets of hyperplanes with a given shape, 
size and orientation. Natural geometrical characteristics of the union of the 
facets, based on Hausdorff measure of the intersections of pairs, triplets, etc., of 
facets form [/-statistics. Building a parametric density from exponential family, 
the limitations for the space of parameters have to be given. In the paper 
basic properties of facet processes are investigated, based on the fact that the 
densities come from exponential family as studied in [5] for the disc process 
with an analogous density type. It is shown that the given class of processes is 
characterized by a repulsive behavior. Conditions for local stability are given 
and a standard simulation algorithm based on birth-death algorithm in Markov 
chain Monte Carlo is mentioned, which enables to observe realizations of the 
facet process. 

As an application of the moment formulas we are interested in the limit 
behaviour when the intensity of the reference Poisson process tends to infinity. 
Central limit theorems for [/-statistics of Poisson processes were derived based 
on Malliavin calculus and the Stein method in [12]. The results were extended 
to the multivariate case in [8] and they apply to the vector of model character¬ 
istics of the Poisson facet process in our setting. In the paper [14] the authors 
give conditions for functionals of Gibbs point processes under which they are 
asymptotically Gaussian with increasing window. Another related work is US, 
where the central limit theorem for the number of intersections in expanding 
window is derived for the stationary planar segment process satisfying certain 
conditions on absolute regularity coefficients. 

In facet processes having densities we restrict ourselves to the model with 
finitely many orientations corresponding to canonical vectors. Submodels of the 
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facet process are investigated where a detailed analysis of the correlation func¬ 
tion yields different asymptotics. When the order of the submodel is not greater 
than the order of the observed [/-statistic then asymptoticaly the mean value 
of the [/-statistic vanishes. This leads to a degeneracy in the sense that some 
orientations are missing. On the other hand when the order of the submodel is 
greater than the order of the observed [/-statistic then the limit of correlation 
function is finite and nonzero and under selected standardization [/-statistic 
tends almost surely to its non-zero expectation. By changing the standardiza¬ 
tion, however, we achieve a finite non-zero asymptotic variance. Even if these 
results are obtained in a special situation with facets of a fixed shape and size 
related to the window size, it is important that they allow us to understand the 
ongoing problems for a possible further investigation of the complex model. 


2 The background of point processes having a 
density 

Consider a bounded Borel set B C R d with Lebesgue measure \B\ >0 and 
a measurable space (N,A/) of integer-valued finite measures on B. AT is the 
smallest er-algebra which makes the mappings x i —> x(A) measurable for all 
Borel sets A C B and all x g N. A random element having a.s. values in 
(N,A/) is called a finite point process. Integer-valued finite measures can be 
represented by systems of points corresponding to their support. Let a Poisson 
point process 77 on B have finite intensity measure A with no atoms, A(P) > 0, 
and distribution P v on M. For a measurable map F : N —> R it holds |1] 


E[F{r))]=e A(B) — [ ... [ F( Ul ,... ,u n )X n (d(u 1 ,..., u n )), 

n—0 U ■ Jb Jb 


(1) 


where we write A"(d(a;i,..., x n )) instead of A(da:i)... A(da: ra ). Further we con¬ 
sider a finite point process 77 on B given by a density p w.r.t. 77 , i.e. with 
distribution 

dP M (x) = p(x)dP r; (x), x S N, (2) 

where p : N —>- R_|_ is measurable satisfying 


[ p(x)dP^(x) = 1. 

J N 

The consequence of ((2|) is a formula 

EF(p)=E[F(r,)p( V )]. (3) 

Let p be a finite point process with density p satisfying 

p(x) > 0 => p(x) > 0 (4) 
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for all x C x. For the (Papangelou) conditional intensity of /x, see pQ, it holds 


A *(u; x) = P( x U x e N, u £ B, u £ x. 
p(x) 

For p(x) = 0 we put A* (it; x) = 0. For ra > 1 we use analogously 

p(xU{«i,.. .,«„}) 


A*(ui, ...,u„;x) = 


P(x) 


ui ,... ,u n £ B\{x} distinct, the conditional intensity of n-th order of /x, Aq = 1 . 
We observe that A* is symmetric in the variables m,... ,u n . The expectations 
of conditional intensities 


p n (ui,... ,u n ; p) =EA*(ui,...,ix„;/x) = E\p(rj U {iti,..., «„})], (5) 

are called n— th correlation functions of the point process /x, cf. [3]. 

For a functional F, y £ B, one defines the difference operator D y F for a 
point process /x as a random variable 

DyF(p) = F(p + 5 y ) - F(p), 

where S y is a Dirac measure at the point y. Inductively for n > 2 and (y i, ...,?/„) £ 
B n we define a function 


D 


Vl, 


F = D 1 D n ~ 1 F 


where D y = D y , D°F = F. Operator D yi Vn is symmetric in yi,... ,y n and 
symmetric functions T£F on B n are defined as 


ir%.fc) = ED;. 

n £ N, TqF = E F(p), whenever the expectations exist. We write T n F for 
T£F. For the functionals of a Poisson process Theorem 1.1 in [7] says that given 
F, F £ L 2 (P V ) it holds 


- 

E[F( V )F(ri)] = EF(r,)EF( V ) + V - (T n F, T n F) n , (6) 

z ' n! 

n= 1 

where (., .)„ is the scalar product in L 2 (A"). We will use symbol [n] = {1,..., n}, 
n £ N. 

For p £ L 2 (Pr/)^ n £ N, it holds 

T„p(yi,...,y„) = ^2 (- 1 )" _|J| p| J |({yj, j £ J};/x) (7) 

JC[n] 

for A n -almost all (j/i,..., y n ), where | J| is the cardinality of J. 
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3 U-statistics 


A [/-statistic of order k £ N of a finite point process /i is a functional defined 

by 

F M = E f(x x k ), (8) 

(xi ,—,x k )en^ 

where / : B k —> M is a function symmetric w.r.t. to the permutations of its 
variables, / £ Li(X k ). Here [& is the set of fc-tuples of different points of p. 
We say that F is driven by /. In this Section basic results on [/-statistics for 
point processes having densities, obtained in [5j are reviewed for later use. By 
the Slivnyak-Mecke theorem m we have 

EF(r)) = [ f{x i,... ,z fc )A fc (d(xi, 

J B k 

For a U- statistic F £ L 2 (P V ) of order k and density p £ L 2 (P V ) it holds 

E F(p) = f(xi,...,x k )pk{x 1 ,...,x k ;ti)X k (d(xi,...,x k )). (9) 

J B k 

In integrals like © it does not mind that p k is defined only for distinct arguments 
Xj. We are interested in higher-order and mixed moments in the following. 

We can use a short expression of formulas for moments using diagrams and 
partitions, see nn, 0 . Let ]~[fc be the set of all partitions {Ji} of [fc], where 
are disjoint blocks and U Ji = [k]. For k = ki -(-+ k m and blocks 

Ji = {j ■ ki d-+ ki -1 < j < k\ -\ -+ ki}, i = 1,..., m, 

consider the partition 7 t = {Ji, 1 < i < m} and let Hfc be the set 

of all partitions er £ J~[ fc such that | J fl,/'| < 1 for all J £ tt and all J' £ a. Here 
| J| is the cardinality of a block J £ a. For a partition er £ k we define the 
function r ■ —¥ M by replacing all variables of the tensor product 

®y=i fj that belong to the same block of er by a new common variable, |cr| is the 
number of blocks in er. 

Theorem 1. Let m £ N, \\T=i F >. e A 2 (P r) ), p £ L 2 (P V ), where Fi are 
U-statistics of orders ki driven by nonnegative functions fi, respectively, i = 
1,..., m. Then 

= E 

«n kl .. 

xp\ a \(xi,... ,x\ a \;p)X{d(xi,.. .,X| CT |)). 

Similar formulas to (flOl) , but under a stronger assumption of a product form 
of the driving function / of F, were derived in [3] using the Georgii-Nguyen- 
Zessin formula pQ. As an application consider processes p a with densities p a 


[ (®iLifi)*(xi,-..,xi er i)x ( 10 ) 

J B M 


E 


n^o.) 
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w.r.t. rj a , a > 1. Let formula m be applied with A a . The term at the highest 
power of a which comes from aX is called the leading term. The rate of the 
leading term is ki. There remains the dependence on a hidden in p, which 

will be investigated later. 

Theorem 2. For a U-statistic 


F(Pa) = 

of order k on a bounded set B we have the following expression for the leading 
term of centered moments, the rate being mk : E[(.F(/i a ) — E F(p a )) m ] = 

X] ( 7 ) J lk /^(Zl. • ■ • . Xlk)pik(x 1, . . . , Xlk] Pa)\ lk (&(xi, . . . , Xlk)) 


' B k 

Proof: We have 


f(xi 


,x k )p k (x 1 ,... ,x k ; Pa)X k {d(xi,.. .,x k )) 


1—1 


E[(F(p a )-EF(p a )r] = J2 

1 =0 


(-1 ) m - l EF(p a ) l (EF(p a )) m - 1 . 


From © it is (E F(p a )) m ~ l = 

\ m—l 

J " f(x 1 ,.. .,x k )pk{x 1 ,.. • ,x k ; Pa)X k {d(x 1 ,.. .,x k ))J 

and from (usd applied to E[F(// a ) z ] we take the term with highest power of a. 
It corresponds to 


a lk / ... ,xik)pik{xi, ■ ■ ■ ,Xik-,p a )X lk {d{xi,... ,xik))- 

J B lk 

These terms come from a &Y\ =1 k ^Cllt w ( lere TI=i fc k a singleton 
with the only partition containing mk blocks all of cardinality one. □ 

Let B C R d , d G N be as in Section 2, X be a point process of compact sets 
m called grains and we denote x realization of X on B , i.e. the union of all 
compact sets. Consider that there is a probability density [5] 

p(x) = c v exp(> • G(x)), (11) 

of X w.r.t. a given reference Poisson point process p of compact sets. Here 
v = (yi, ..., i/d) is a vector of real parameters, c v is a normalizing constant, 

G(x) = (G 1 (x),...,G l! (x))eI d 
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is a vector of geometrical characteristics of x. In the exponent of m there is the 
inner product in R d . The largest set of v such that exponential family density 
HID is well defined is 

{v G : E[exp(z/ • G(?y))] < oo}. 

For m G N denote 

= in;: .i i x .)./>";. 

the vector of m—th differences. 

Proposition 1. Consider the probability density m- Then for the correspond¬ 
ing conditional intensity of order m G N and a realization x of X it holds 

A^(yi, • • •, y m \ x) = exp(> • Q m G(x)), (12) 

where y 1 ,... , y m G B \ x are distinct, 

Q m G(x) = D™ i ... i „ ra G(x) 

+ E + £ D Vi G(x). 

G{l,...,m} l<i<m 

The intensity of the reference Poisson process depends on a specific model, 
see [9] for interacting discs. 

4 Facet processes with density 

Here we consider processes of interacting facets in R d , d G N, with their natural 
/7-statistics. Let 

Y = Bx (0,6] xS d_1 , (13) 

where b > 0 is a size parameter, S^” 1 is the hemisphere of axial orientations in 
R d . For a point y G Y, y = (z,r,c/>) represents a facet which is a subset of a 
hyperplane A = A(z, (f>) through point z having normal orientation </>. Then 

y = {s G A; dist{s , z ) < r} (14) 

for a distance dist in A. The Poisson process y on Y has intensity measure A, 

A(d(z, r, </>)) = x(-)d^G(dr)u(</>)(d^>), (15) 

where Q is the size distribution, a probability measure on (0, b], v is a probability 
density w.r.t. Lebesgue measure on § d-1 , y is a bounded intensity function of 
facet centres on B. 

Let k G {2,... ,d}, we consider the intersection of any k facets such that 
the corresponding hyperplanes are in general position, cf. (13], p.133. The X k - 
measure of all such fc-tuples is equal to A (Y) k for each k since we have a density 
v in m- For the Hausdorff measure HP in R d of order j and x G N (on Y) we 
put 
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(16) 


G 1 (x) = ^H d - 1 ( y ), 


g 2 (x) = 1 E w d - 2 ( yi n y2 ), 


(y i,»2)ex^ 





E ^(nti 1 ^), 


G d (x) = i E H °(nf = i yi ). 

(lit 'll j)Cy^ 


(»,-,w)£xJ 


We deal with 


G(x) = (G 1 (x),...,G d (x)), 


(17) 


here Gj is [/-statistic of the j —th order, j = 1,..., d. The facet process p is 
defined by the density CD 


p(x) = c v exp (v ■ G(x)) 


with respect to p, where v = (v \ 1 ..., z^) £ R d . For x = {ui,..., u n } we have 

Gj(x) < const.n^b d ~ :i , 

which together with the assumption i/j < 0, j = 2,... ,d proves using CD that 

P £ Li(P v ) n L 2 (P V ). 

In fact the assumption that the orientation distribution V ( d(f >) has a density u 
can be weakened to the assumption that V has at least d atoms, which we will 
use later. 

Example 1. Specially for d = 3 the facet process may serve as a model for 
platelike particles in materials micro structure of metals. Here G\ yields the 
total area of all plates, G 2 is the total length of intersection segments of pairs of 
particles, G 3 is the total number of intersections of triplets of particles. The size 
of negative parameter v 2 or 1/3 gives the measure of neglection of intersections 
(repulsion) of the corresponding type. If v is a zero vector, then the facet process 
is Poisson (with no repulsion). 

We will consider special types of facet processes p^ k \ k = 1,..., d, such that 
in CD we have Vk < 0 while Vj = 0, j ^ k. We say that ^ k ) i s a submodel of 
order k , especially p^ is a Poisson process and p^ is a Poisson process if and 
only if Vk = 0, k = 2,..., d. 
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Local stability of a facet process is defined by the existence of a > 0 such 
that 

A* (it; x) < a, x S N, u £Y \ {x} (18) 

and under given conditions the processes pS k \ k = 1 ,,d are locally stable. 
Simulation of a facet process p is available using the birth-death Metropolis- 
Hastings algorithm [5]. 

In |9j the process with conditional intensity A* is called attractive if 

A*(y;xi) > A*(y; x 2 ), 

for any Xi,x 2 G N, X 2 C xi. It is repulsive if 

A*(y;xi) < A*(y; x 2 ), 

for any xi,x 2 £ N, x 2 C xi. In the case of sharp inequality we say that the 
process is strictly attractive, strictly repulsive, respectively. 

Proposition 2. The facet process /Z 1 ) is neither strictly attractive nor strictly 
repulsive for any u\ £ R. The facet processes p^ are repulsive for k = 2,... ,d. 

5 Asymptotics with increasing intensity 

5.1 The Poisson case 

Generally on B C R d , A as in Section 2, for l > 1 and i = 1,..., l let Aq £ N, 
/W £ Li(A ki ) be symmetric functions. Consider Poisson processes r] a with 
intensity measures A a = aA, a > 1. Following [8] P-statistics 

F^(Va)= E / (<) (* X ki ) (19) 

(xi,...,X k .)eVa^ 

are transformed to 

F« = a -( fc -5)(FW - EfW). (20) 

The asymptotic covariances are 

Cij= lim cov{fW,fW)= f T 1 F^(x)T 1 F^(x)X(dx), i,j£{l,...,l}. 

a ~> °° JB 

(21) 

The convergence under the distance between Z-dimensional random vectors P, Z 

d 3 (U,Z) = BUT>\Eg(U)-Eg(Z)\, 
g&n 

where Tl is the system of functions h £ (7 3 (R ( ) with 

, d 2 h(x) , d 3 h(x) 

max sup ---— <1, max sup -----— <1 

l<ii<i2<l xGR i ' dx^dx^ ' l<ii<i 2 <*3<i xG R( dx^dx^dx^ 
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implies convergence in distribution. Based on the multi-dimensional Malliavin- 
Stein inequality for the distance d 3 of a random vector from a centered Gaussian 
random vector Z with covariance matrix G = [5] show that under 

the assumption 

[ |TiF«| 3 dA < 00 , i = (22) 

J B 

there exists a constant c such that 

d 3 ((F^,...,F^),Z)<ca~K a > 1. (23) 

Example 2. For the Poisson facet processes rj a , a > 1 on Y \13\) with intensity 
measure aX fW and the U-statistics Gj(r]),j = 1,..., d, in it 161) we obtain that 

TiGj(x) = yyy J ■ ■ ■ na;)A J ” J (d(yi, 1 )). 

The finiteness of the intensity measure A in m and the boundedness of the 
facets guarantee that all integrals it jj|) and \21\) are finite. Thus for the random 
vector ( Gj(p a ),j = 1 ,d) both the central limit theorem when a —> 00 and the 
Berry-Esseen type inequality fUj) hold. 


5.2 The non-Poisson case 

Let facet processes p a , a > 1, have densities 


Pa(x) = c Vt a exp (u • G(x)) 


(24) 


w.r.t. Poisson processes rj a with intensities A a = aA, respectively. Here G'(x) is 
given in m and uj < 0, j = 2,..., d. We investigate [/-statistics and submodels 
of the same order d — k, from formula ([9]) we have 

E G d . k (^~ k) ) = (25) 


~.d—k 


0 d-k)\ 


k = 0,..., d — 2. 


f H fe (nt^)p d _ fc (yi,„ 

! Jyd-k 


Lemma 1. For the processes pi k> the correlation function pu-k , k = 0,.... 
d — 2, has form 


Pd—k{yii • • • > Ud—kj Pa — 


A H 

B(aY 


(26) 


where A(a) = 

00 „n f 

E 


! J Y '"Jy CXP \(d~k)\ 


Vd-k 


n —0 


E 


uk (ntiXi) 


{xi,...,x d - k )C{u 1 ,...,u n ,y 1 ,...,y d -k} 
xA"(d(ui,.. 
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B{a) = Y 


n —0 


n/- 

• M 

n\ J Y 

Jy \ 


k'd-k 

(d - k)\ 


E 




(xi,...,x d -k)C{u 1 ,...,u n } 


xA"(d(ui,.. .,u n )). 

Proof: Using formulas (JSJ,(]TJ) we obtain 


Pd-fe(£/i, - - ■, 2/d-fc; fe) ) = c„ iQ e aA(Y) ^ 


n=0 


a'" 

n! 


/y 


exp 


Vd-k 

(d-k)l 


E 


I A"(d(u 1 ,... ,«„)). 


(xi,...,x d -k)C{u 1 ,...,u n ,y 1 ,...,y d -k} 

The normalizing constant can be expressed from 


1 = Ep(r] a ) = c„, 0 e aMY) Y 


n =0 


a 

n\ 


exp 


Vd—k 


(d-k)\ 


Y H fc (n£i**<) | A"(d( Ul ,...,« n )) 

(x 1 ,...,x d _ k )c{u 1 ,...,u n } 

and the result follows. 

We obtain asymptotic results in a special model. Let 


Y = [0,6]“ x {26} x {a, i = l,...,d}, 
where e, are canonical unit vectors. A facet is a set 


□ 


(27) 


y = (z, cj)) = {s G D; max \z t - s*| < b}, 

i 

z denotes the centre and D is hyperplane with normal orientation <fr. That means 
facets have the same fixed size and shape and any non-parallel facets intersect. 
In the case of d — k facets with different orientations we have bounds 


b k < H fc (n fl k Xi) < (2 b) k . 


(28) 


In the intensity A (usd we have Q = 52b, the orientation distribution V is uniform 
on {e*, i = 1,..., d}. Then for u* = (zi, </>,;), i = 1,..., n, it holds 

A"(d(u 1 ,...,u n ))= (29) 


x{z\)Az \... x(z n )dz n —Y ■ 


E n\ 

- I 

niH-h n d =n n\\ . . . Tlrf)- 


^^(d^i,... ,d<t>n)- 


Denote T = / [ob]d x(z)dz. 
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Remark 1. In the special model it can be shown that for any l € {2, ...,d} 
if Vi = 0,* l and vi > 0, then p ef L 1 ^^), thus the conditions applied to 
parameters are not only sufficient, but also neccesarry conditions for density 
existence. 

Lemma 2. Denote as in Lemma Q] formula 12(A ) 

Pd-k(yi, ■ ■ ■ ,yd-k', ^a~ k) ) = 4t-t> k = 0,...,d-2. 

B{a ) 

Then we have for yi, ..., y d -k with different orientations 


B(a ) > ^ 


n =0 


aTV 


E-E 


niH-h n d =n n\\ . . . Tlrf! 


x exp 


d—k 

V d - k (2b) k Y n n h I > 

d-k i —1 


A(a) < Y 


n —0 


x exp 


oTV 

~d) 


( 


E-E 

/ 


1 


iH-hn d =nn 1 ! . . . rid! 


Vd-kb 


\ 


d—k 


e n^ 


V d—fc _ rm *=i 


c[. 

where d is the number of (d — k)-tuples (with different orientations) 


{xi, ..., x d -k} c {ui, ...,Un,yi,..., yd-k} 


such that there is at least one of yj among xi,.. . ,x d -k- 

Proof: It follows from Lemma U inequality (1251) and integration w.r.t. (12U1) . 
Here rij are interpreted as the number of ui, l = 1, ... n with orientation e^, j = 
1,... d. □ 

Lemma 3. Under the notation from Lemma [H it holds 

B(a ) > exp(— aT(d — k — 1)), k = 0,..., d — 2. 
d 
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Proof: Substituting /? = we have 


S' 3 " E 

n—0 ni~\-...n d =n 


oo oo 


d—k 


n\\. ..n d \ 


exp(^ d _ fc (26) fc e ik>= 

{li}*l?c[d\ j =1 

d—k 


^ w _ /oni+...+n d _ __ 

e • • ■ e ni! ,.. nd! ex p( z/ rf-fc(2&) ^ n%)= 


m=o nd—O 


Mlzfcld] j=i 


E ••• E tn—I—r ex P 


/ 


m=0 nrf_i=0 


n\\. ..n d - 1 ! 


d—k— 1 


/9exp(i/ d _ fc (26) fe e n %) i >< 

{^}?=o fc_1 C[d-l] J =1 


/ 


x exp 


d—k 


v d -k(2b) k ^ IT n b I — 

(b'}?=i fc C[d-l] i=l 


d—k 


EE EE / Qni + ...+n d _ 1 

E--- E ni \...n d . 1 \ eMVd - k{2b) E IW> 


«i=0 ra d _i=0 


{b}ti k c[d-i]J =1 


d—k 


EE EE /Qni + ...+n d _ fc 

E - E »,!,..n j _ t ! exp( ‘ / -‘- t(2l,) '' E n%) = 


ni=0 rid—k—0 


{^}?=i fc C[d-fc] J =1 


d—k 


EE EE /Qni + ...+n d _ fc 

E-" E Ej-—- T exp(^_ fe ( 26 ) fc n n ‘) = 


m=0 nd-k= 0 

oo °° on, -I- ,.+n d _ fc _ 1 


E- E 

rai=0 ra d _ fc _i=0 
OO 


n\\. ..n d - k \ 

pn 1 + 

ni\ ■ ..n d -k- 1 ! 


i= 1 

\ k I— rd — k — 


exp^-^nt*- 1 ^ 


> 


E- E 

rai=0 ra d _fe_i=0 


/3 ni+ - 


, + Ud-k-l 


nil... n d -k- 1 ! 


— e /3(rf—fc—1) 


□ 


Lemma 4. There exists a constant R < 0 independent of a such that 

A ( a ) < Ra 

B(a) ~ 


Proof: Substituting j3 = we estimate A(a) 


(30) 


A(a)<y2/3 n 'V] —j--exp (v d _ k b k C d -k,d), (31) 

Eo n 1 + .E d =n niL '- nd! 
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where 


l i 

Cl,m = } ' 'y ' n jk + Ij 1 ^ — m - 

i=1 Oi>---J<}C[to] fc= 1 

^l,m = C'l—i m—l'ft'm H” ^Z,m—1? 2 ^ 771, 

Cm,m — ^ra—l,ra—1 (1 H - ^m) • 

Now we decrease the number of sums in (f3TT) : 


EE EE /Qni + ...+n d 
V] • • • V] —;-r exp {v d - k b k C d _ kA ) = 

Eo ^o n i ! "- nd! 

^ni + ...+n d 

^ ^ ^ ^ “ j “ j~ ^Xp(z/ £ f_^& (Cd—fc— l,cZ_l Tld “f Cd—k,d— l)) — 


ni\...n d \ 

m=0 rid—O 

00 00 oni + ...+nrf_i 


E E 

" 1=0 n d _ i=0 


nj.. .n d _i! 


oo oo 


yjni + ...+n d _i 


ni!. ..n d - 1 ! 


E E 

"1=0 "d-l=0 

p p yJ"i + ...+"d- 


exp(/3e Vd - kbkCd - k - 1 ’ d - 1 + v d - k b k C d -k, d -i 
exp {/3e Vd - kbkp + u d _ k b k C d _ kid - 1 )+ 


E- E 


m!... nd-i! 


exp(/3i 


d-kb k Cd-k-l,d-l 


+ v d -kb k C d -k , d -1 


"1=0 "d-i=0 

p p pn 1 + ...+n d _ 1 


E E 

"1=0 "d_i=0 


ni!.. .n d -i! 


exp(/3e iyd_fchfcp + v d - k #C d - k ^{). 


)< 


)- 


In the last step (p+l) d_1 terms of infinite series were changed. Then we estimate 
the exponential part of series from above by exp(/3e Vd ~ kP ). 

Each term of both finite series tends to zero when divided by e^( d_fe_1 ) (as /3 
tends to infinity), so that the whole series tends to zero. This step is repeated 
till we get 


E- E 


on 1 + ...+n d _ k 


P , i I 

^ n ni!.. .n d _ fc ! 

" 1=0 n d _ k —0 


OLr 


>E- E 

ni=0 n d _ k = 0 
OO oo 

• E- E 

" 1=0 n d _ k _ i=0 


0 n 1 + ...+n d - k 

nd . ..n d - k \ 


exp {v d _ k b k C d _ k ' d _ k ) < 

—k—1 d— 

exp | v d - k b k ^2 n i n d-k + E 


—k 


Vd-kb \ni+...+n d _fc_i 


(Pe Vd ~ kb ) 


i =1 


77i! . . .rid-fc-i! 


■ exp(/3e I/d_t:b ^ i=1 n * +1 )) < 
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a t 


°° 00 (a p Vd-kb k 'ini + ...+ri < 1 _fc-i , 

E- E (ge ni ,.. ) e*p(^-^) + 


ni=0 rid—k—l—O 
9 9 


._. ._. (RpUd-kb k \n 1 +...+n d _ k _i , , 

E ^ , -i- eMPe^^- n ^y 

n\=0 rid -EUo TT-i! . . . TT-d-fe-l! 

“pE". E -—-1-T- exp(fie Vd ~ kb q ), 


n i =0 n d - k - 1=0 


nfi ... .n d _ fc _i! 


where a p = exp(/3e l,d ~ kb p k). Again we changed (q + i)**-* -1 terms of infinite 
series and estimated exponential part of the series by exp(/3e Vd ~ kg ). Both finite 
series tend to zero when divided by ePG-b- 1 ) and the infinite series is equal to 


exp(/3 e Vd ~ kbkp k + pe Vd ~ kbk (d-k- 1) + fie Vd ~ kbkq ). 


We have to select p and q so that 

Rl = e ^d-kb k p k + e v d -kb k (rf _ k _ i) + e "d-kb k q _ ( d _ k - 1) < 0, 

which is possible for any u d -k < 0. Then we reverse the substitution /3 = s j- 
and the ratio tends to zero for a —> 00 . Its convergence is not slower than 
e Ra where R = . □ 

5 a 

Theorem 3. For a —> 00 it holds 

EG d - k (pi d ~V) ->0, k = 0,... ,d — 2. (32) 

Proof: In (E5l) it suffices to show that for yi,... ,y d -k all with distinct orien¬ 
tation we have 


lim a d k p d -k(yi, ■ ■ ■ ,y d -k\ l^a fe) )=0. (33) 

a—> 00 

This is the consequence of Lemmas [3] and 0] Then by the Lebesgue dominance 
theorem formula (13^1) follows. □ 

Remark 2. Since the functionals G d -k are nonnegative Theorem^ says that 
asymptotically the processes t L a~ k are degenerate (for any v d -k < 0 ) in the 
sense that there are no intersections of (d— k)-tuples, i.e. some orientations are 
missing. For k = d — 2 e.g. all facets tend to he parallel (with any orientation). 
Additionally when the assertion of Theorem [3 is valid for p~ d ^ then it holds 
also for expectations E Gj(p^a k ^) with j = d — k + 1 ,... ,d. 

Remark 3. The behavior of Metropolis-Hastings chain for simulation of a real¬ 
ization of p d ~ k with fixed parameter a depends on v d -k in the model with finitely 
many orientations. For \ v d -k\ large typically it converges quickly to a realization 
with some missing orientations. However when v d -k is close to zero it does not 
converge quickly at all. This property does not take place when the orientation 
distribution is absolutely continuous w.r.t. spherical Lebesgue measure. 
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6 [/-statistics of the order smaller than the sub¬ 
model 

In this section we continue to study the model from previous subsection, see m- 
Asymptotic moments (when a —> oo) of functionals Gd-fc(x) in the submodel 
n\-f \ where k = 1,..., d — 1; d > 3, will be investigated. Here the vanishing 
property © is not expected, see the following Table where crosses mean that 
the expected value is non-zero. 


a —>■ oo 

H-statistics Gj 

Submodel 


EG d _r 


eg 2 

EGi 

u (2) 

0 

0 


0 

X 

U (3) 

pa 

0 

0 


X 

X 







U (d_1) 

pa 

0 

0 


X 

X 

U (d) 

pa 

0 

X 


X 

X 


As in the previous Section we will need first some limits of correlation func¬ 
tions. 

6.1 The limit of correlation functions 

Proposition 3. a) For any yi,... ,y d -k G Y with different orientations we 
have 

k 

lirn p d -k{yi,- ..,y d -k;ijl d) ) = -r, k = 1. (34) 

a—too CL 

b) For all sets of2(d—k) arguments which fullfil that each facet of {z/i,..., yd-k} 
has different orientation and the same applies to {y d -k+ i, •.., 2 / 2 (d—fc)} a nd there 
is l common orientations within these two groups, then 

Pw-k)(yi, ■■■, V2 vi d) ) = d +l ' 

c) Consider 2(d — k) — 1 arguments such that each facet of {y \,..., y d -k} has 

different orientation and the same applies to {y d -k+ i, • • ■, 2/2(d—fe)—l? Vi}- Let us 
have l common orientations ei ,..., within these two groups of arguments 

not including orientation ei of y\, then 

.. , (h)\ 2 k — d -f-1 — 1 . . 

lim p 2 (d-k)-l{yi, ■ ■ ■ ,y 2 (d-k)-l,ffa ) = --7-■ (36) 
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Proof a) Consider fixed facets yi, ■ ■ ■ ,yd-k each associated with a different 
normal orientation vector, without loss of generality e\,..., e d -k- Then the cor¬ 
relation function can be expressed analogously to Lemma |T] as follows: 


Pd-k(yi, ■ ■ ■, Vd-k\P { a ) ) 


Ma) 

B( a y 


A(a) = 



A n (d(ui,.. 


A”(d(ui,.. .,u n )). 


It holds A(a)/B(a) = A(a) /B(a). where using (HHil) we write A(a) = e Ar( d Lx 


£(% 


n —0 


E 


i 


n 1 + ...+n d =n 


ni! ■ ■ -n d \ 


exp | v d e n ni ) , (37) 

Dc[d-k] l£[d]\D 


= 1) V — 


B(a) = e d 


n—0 


aT 

d 


n 1 + ...+nd=n 


ni\...n d \ 


exp 



From Lemmas [5] and [6] below, (l34l) follows. 

b) We need to compute P 2 ( d -k) (yi , • ■ •, V 2 (d-k)\ PaP) for all sets of arguments 
which fullfil that each facet of {y i,... ,y d -k} has different orientation and the 
same applies to {y d -k+ 1 , ■ • •, 3 / 2 (d-k)}- Without loss of generality we need to 
consider only situations where facets y 1 ,..., y d -k have orientations ei,..., e d -k 
and facets y d -k+ 1 , • ■ •, 2 / 2 (d-k) have orientations ei,..., e u e d -k+ 1 , • ■ • e 2 d - 2 k-i, 
where l = max(d — 2 k, 0 ),..., d — k, probability of both sets of facets having l 
common orientations can be expressed by multinomial coefficient 


(d — k )\ 2 / d 

d 2 fo- fc ) \l,d —k — l,d —k — l, 2 k — d + l 


(38) 


It holds 


—a(d— 1 ) 


hm 2> 2 (rf—fe) ( 2/1 , - - •, 2 / 2 (d—fe) 5 Ai£ d) ) = hm ---V 

a—¥ 00 v ' v ' a—¥00 rl z ' 


a->oo d 2 -—< Tl\ 

n =0 


exp 


d\ 


E H 0 (n? =1 *0 | A"(d(in,...,«„)), 
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where we get -—^—- from calculations of the limit B{a) in Lemma U For 
selected l the limit of the correlation function can be further expressed in the 
form 


-a(d-l) 

lim ---- 

a—too d 


E-E 

m=o nd=0 


ni! ...n d \ 


x exp 


^ e n ’ 

D(Z[2d—2k—l] m^[d]\D 


2 I-Dn{i,...,z}| 


(39) 


The expression in the exponent can be bounded from both sides 


Vd 2 Z ^2 n m < 

D(Z[2d—2k—l\ m£[d]\D 

^ e n 2 |jDn{i . i}| < 

Dd[2d—2k—l] m€[d]\D 

Vd y. n m . 

D(Z[2d—2k—l] m(E[d]\D 


For the both bounding series the limit (15^1) is equal to 2fc ~ rf d+i since the series 
are in the form as A(a) in Lemma [6] 

c) In the previous case it can be seen the numerator is equal to number of 
unused orientations among facets y\, ■ ■ ■ , V 2 (d-k) and the same applies to this 
calculation except that there is one extra common orientation of mi. □ 

Remark 4. In b) the limit is equal to zero for l = d — 2k. In this case all 
orientations up to the order d of the submodel are exhausted. This corresponds 
to the situation from Theorem [3 where in 1,9.9)) also all orientations up to the 
order d — k of the submodel are exhausted by y i,... ,y d -k o,nd the correlation 
function tends to zero. In the opposite case (in Theorem^ when l > d — 2k is 
admissible) the limit of correlation function is nonzero. 

In the following two lemmas we write for simplicity v instead of v d . 

Lemma 5. It holds lim^oo B{a) = d. 

Proof: We want to examine series in form 

~ ~ a ni+...+n d 

/■■■/, —1 - r exp^ni .. ,n d — a{d — 1)) = (40) 


E E 

rii=0 nj_i= 0 




n\\.. .nd-T- 


exp(cte 


vn 1 ...n d _ 1 


- a{d — 1)), 


where v < 0 and d > 3. This form fits the B{a) in (l37l) after substituting a for 
pf. We divide indices in the sums into three subsets 

D\ = {ni > \fa A ... A rid_i > y/a}, (41) 
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D 2 = {m = 0 V ... V rid-i = 0}, 

D 3 = {n\ < y/a V ... V rid -1 < y/a} \ D 2 . 

Then we use one of Chernoff’s bounds for Poisson distribution, which says that 
it holds 


‘ - ^- s (esY 


7T e_s — e_ 

' k\ 


k =0 


P 


,t < S. 


Firstly, we sum over D 2 . Let m = 0, then from (140]) we have 

1 n 2 + ...+n d _ 1 


E E 


n 2 \. ..n d - 1 ! 


exp (a - a(d - 1)) - e a-a(d-l)+a(d- 2 ) = L 


n 2 =0 —1=0 

We continue by induction with m 0,..., n ^-1 y^ 0, n^. = 0 where we get 

ni + ...+nfc_i+nfc + i + ...+nd_i 


OO OO 


EE E E 

ni = l n fc _i=l n fc+ i=0 n d _i=0 


a" 


m!...nfe-i!nfe + i!.. ,n d - 1 ! 


exp(a — a {d — 1)) = 


= e _° (fc _i)( eQ _ 1 ) fc _i ^ L 


By calculating all d — 1 options we explored all combinations of indices in D 2 
and we conclude that the sum over D 2 tends to d — 1. 

Secondly, we sum over £> 3 . 


E 

{n 1 ,...,n d _ 1 }eD 3 

E 

{ni,...,n d -i}£D 3 


a ni+...+n d _i 

ni!...n d _i! 


exp(c 


vn\...rid -1 


— a(d — 1)) < 


m + ...+n d - 1 


n\\.. .rid- 1 ! 


■ exp(ae y — a(d — 1)) < 


e ae “ E 

{ni,...,n J _ 1 }eD 2 U 03 

/ L%/“J 


ni+...+n d _i 


n\\... rid- 1 ! 


exp(— a(d — 1)) < 


“ME 

\m=0 


mi 


E 


nd-V 


< 


e ae “{d- l)e 


n d -x=0 
-a M L ^ J 

(L^j)WHj’ 


where we used the principle of inclusion and exclusion (for probabilities of Pois¬ 
son distribution) and the Chernoff’s bound. Then we examine logarithm of the 
previous term 

a{e v - 1) + log(d - 1) + [yfa\ log(ea) - [v^J log(LVaJ)- 
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which tends to — oo and so the sum over D 3 tends to 0 , when a —> 00 . 

Finally, we examine the sum over indices in D\ in order to show that it tends 
to 1. Let e > 0 be arbitrarily chosen. We choose 71 , so that 

| exp(ae I '^' v/ “^ ) — 1| < e,Va > 71 . 


In the next step we choose 72 , which fullfills the following condition: 


ea 




(L^J)Lv^J 

Then we estimate the series from both sides 


< e, Va > 72 . 


E E 

ni = |V“1 nd-i = \y/a\ 


rii + ...+ra d _i 

“_ e -“(d- 1 ) 

ni\. ..n d -\\ 


< 


(42) 


E E 

ni=f\/al n d _i = [v/al 


ni+...+n d _i 

__ c aG 1 d ~ 1 —a(d— 1) ^ 

nx\. ..rid- 1 ! 


E E 

™i=T\/“l ri d - i = |V“l 


a 


ni + ...+n d _i 


m! • ■ • nd_i! 


e -a(d- l)( 1+e ) < (1+e). 


Investigating the lower bound it can be seen that 


OO OO 

E - E . 

n 1 — \ y /a] n d _i = [ v /a] 

1- E 

{n 1 ,...n d _ 1 }GD 2 UD 3 


a ni + ...+n d _i 

nil... rid- 1 ! 
ni!...n d _i! 


e -a(d- 1) _ 
, e -a(d-l) > 


Lv^J 


l-(d-l) E 


> 1 - (d- 1 ) 


n =0 


-a_(fEEL 

(L^j)L^j 


> 1 — (d — l)e, 


where we used again the principle of inclusion and exclusion and the Chernoff’s 
bound. Thus the sum over D\ can be enclosed by bounds which are arbitrarily 
close to 1. We conclude that the overall sum CU1) tends to d. □ 

Lemma 6. It holds lim^oo A(a) = k. 

Proof: We examine the expression for A(a) in (1371) substituting ^ = a, in the 
form 


a ni+...+n d 

V ■ ■ - E —i-T exp {vQd-k.d -a(d- 1 )), 1 < k < d - 1 , (43) 

„ ^ n 1 \...n d \ 

n i=0 rid —0 

where 

Q s ,t = e n »«,o< a < t , 

FC[«] Je[t]\F 
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n ie0 n ; = 1- It holds 


Qa,t ^ L 

Qs,t = Qs,s^s+1 • • • rif, S < t, 

Qt,t = Qt-i,t-i(nt + 1 ), 


and it follows 

~ “ a ni+*-+n d 

V • ■ • V —j- r exp {uQd-k,d - a(d - 1)) = 

n i =0 rid —0 


= E E 

ni=0 n d _ 1 =0 


a n 1 + ...+7i d _ 1 

nil. ..n d - 1 ! 


exp(c 


d—k,d — l n (d 


(d- 1)). 


(44) 


For the same subsets of indices as in dSD we use the fact that current series are 
bounded from above by the corresponding ones. Therefore the arguments from 
Lemma [5] remain the same for D 3 , and also for D\ since (l44l) is ordered between 
the first and second expression in (1421) . 

The subset D 2 has to be considered where we proceed analogously to the 
proof of Lemma 7, but in reverse order from n d ~i to n\. The fact that in all 
terms of Qd-k,d- 1 there are factors n d -i, ■ ■ ■ nd-k+i is used. Setting rid-t to 
zero we get 


OO OO OO OO 

E- E E E 


1 +...+nd _ t _ 1 + 71 ^ _ t -|_ 1 +...+rid - 1 

raj.!... n d -t-i } -nd-t+i } - ■ ■ • n d - 1 ! 


711=0 Tid — 1 — 1=0 Tl^ — £+1=1 71 — 1=1 

x exp (a - a(d - 1)) = (e a - ^-^^-^ 1 )^- 2 ) 1( k > t> 


OOOOOO OO 

E- E E E 


a 


, n 1 +... -\-rid -1 -1 -\~rid _ t +1 + • • •+ n d -1 


ni!... n rf _t_i!n d _t + i!... n d - 1 ! 


711=0 7ld — t — 1=0 71d — £+1=1 7ld — 1=1 

x exp {ae v - a{d - 1)) = (e a - iy-i e a(d-t-i) e av-a(d-i) 0] k<t, 


where the second expression is an upper bound. Thus the sum over D 2 tends 
to k — 1 and we conclude that the overall sum (H51) tends to k. □ 


6.2 The asymptotics of moments 

In this subsection first and second moments of functionals G d -k{yf) in the sub¬ 
model [in ' 1 , where k = l,...,d— 1, d > 3, are studied. We will need the 
following Lemma. 

Lemma 7. Let two sums of finite length XEi an d fulfill 

K L 

lim V nf ‘= lim y' m\ a \ 

a—100 1 a-t 00 1 

i=l i=l 
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(a) _ c ),i (a) r / ? 


7^’ 




lim cfl = Cj,! G (0,oo), lim c-^ = Cj i2 G (0, oo), 

a—t oo ’ a —>00 ’ 


lim d,| a L ' ) = d,;,i G (0,oo), lim d)°o = d i]2 G (0,oo). 

a—>-oo ’ a—>-oo ’ 

Lei t/iere exist ao > 0 , 7 < 0 and £ > 0 such that for all a > ao and all i it 
holds 


( a ) 


t/ien /or p > 0 it holds 


(a) 1 

q, 1 - ci,i| 

< Ce 7 °, 

lc (a) 

l C 2,2 

4“i - di,i| 

< Ce 7Q , 

i45 


lim a p 

a—too 


K „(“) i j(“) 

y^ fiii _ y^ 4i 

Z-z la) ,(a) 

z—1 H ,2 z=l u i,2 


= 0 . 


Proof: It holds 


lim a p 

a—too 


^ J“) 

L 

j(«) 


K So.) 

i j(“) 

c z,i 

Z- ( 0 ) 

z= 1 c z,2 

-E 

2=1 

a 2,l 

< lim a p 

a—too 

y— c z,i 

Y^ 4,i 

d (o) 

a i,2 

r-f Ci , 2 - eT« 
2=1 ’ 

*-? di, 2 + 

2=1 ’ 


yd linia-^oo a p 4“/ y-^ lima^oo 


Ci 2 


3i,2 


< 


if 

E 


lima-^oo a p Ce 7 “ 


c; 2 


E 


lim a ^oo a p fe ia 


M,2 


k l , 

Q,1 y^ dj,l 

a->oo c Z_z ^ 

2=1 ’ 2=1 ’ 


lim 


= 0 . 

□ 


We denote 

4= [ 


ITfc (nf = i fe (s i ,26,e i ))x(si) •.-x(sd-fc)dsi...ds d _ fe , (45) 


/([0,b] d ) d - fc 


4 = / H fc (nt 1 fc (a < ,26 ) e i ))H fc (nJ!r a fc (s i+ d- fc -i,26 ) e i ) n (si,2&, ei ))x 

d ([ 0 , b ] d ) 2 ( d - fe )- 1 

xx(si)dsi,... ,x(s 2 (d-fe)-i)ds 2 (d_fe)-i, 

for facets yi = (s*, 2&, e*), i = 1,... ,d — k, with different orientations. Because 
of symmetry 4; do not depend on the choice of these orientations. From 
((Ml) it is 

0 < T d ~ k b k < 4 < T d ~ k (2b) k . (46) 
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Theorem 4. It holds for k = l,... ,d — 1 : 


lim 

a—too 


G d - k (p {d) ) 

n d—k 


h 

d d ~k 




(47) 


lim — 

a—too a 


1 

2 (d-k)- 


1 


var 


G d -kU d) ) 


, (d- l)(d- k) 2 

±k ^ 2 (d-fc)-l 


/ d -2 
U - fc - 



Proof: First we show that 


(48) 


Hm E G d . k U d) ) = I k /d-l\ 

a->oo a d ~ k d d ~ k \d — k) 

From (0 and m we have 

EG d . k (p^) = 

= j^y j M. k (n d = 1 k y l ) Pd -k(yi, ..., y d - fc ; p^)X d ~ k (d( yi , ■ 


(49) 

(50) 

• i Vd-k))- 


Finally in order to obtain (l49l) we use (J34]), the fact that d< ' d ~ 1 ' > d d- d ~ k+1 ' > is the 
probability of facets having different orientations and the Lebesgue dominance 
theorem in (l50l) . 

Further we evaluate the second moment of the tl-statistics according to 
Theorem 1 using notation 


m k ( yi ,...,y d _ k )=m k (n d i k yi ). 

By — we stress that only terms with higher or equal rate of convergence than 
a 2 (d—fc)—l are expressed: 


G 2 d _ k { { u^) = Y, 


i.M 


((d-k)l) 


(51) 


)Y M 


a ^Tld.-k,d-k 

(yi, ■ ■ ■ > yw\)pw\ (yi, ■ ■ • > vw \; p {d) )Kd{y\, ■ ■ ■ > vw\)) 


— J\ + J 2 , 


where 


Ji = 


a 


2 (d-k) 


[‘(rtiV)H‘(n&ii y*)x 


((d — /c)!) 2 J Y 2(d-k) 

xp2( d ~k){yi, • • •, y2(d-kf P (d) )>? {d ~ k) (d(yi, ■ • ■ >?/2(d-fc)))> 

T Q ^ (d k)‘ f TTir k/yxL—k \|n]-fe/p^ 2 (d—fc) —1 „ \ 

J 2 - ((d _ fc)!) 2 Jy2(d — fc) — 1 H ( i=1 yi)H 2/1 0 yi)X 

xp 2 (d-k)-l(yi,- ■ ■ , y 2 (d-k)-l: P { a d) )^ 2id ~ k) (d( yil . . . , 2 / 2 (d-fc)-l))' 
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Here (d — k) 2 is the number of possible selections of a common element in the 
both functions. The first term of (l5lT) tends for a —> oo to 

II (d- 1\ 2 

d 2 (d-k) y d — kj 

From the Lebesgue dominance theorem in J\, using (1551) . (1551) . (1551) : 

EG d _ k (/di -*) _ dfc y- / d \ 2 k — d + l 

a->oo a 2(d-fc) d 2 G~ k ) — k — l,d — k — l,2k ^ d + l) d 

II y^Y d ~ 1 ) Ik f d ~ 1 \ > 

~ d 2 (d-k) 2_, \j )( i-k-l,d-k-l,2k-d + l- l) ~ d 2 G~ k ) \d-kJ ' 


Here the last equation holds because both sides represent (up to a multiplicative 
constant) number of possibilities how to choose two sets of distinct d—k elements 
from total d — 1 elements. This term cancels out with the squared expectation 
limit in (l49l) . therefore 


lim 

a—too 


var Gd-kj^a' 1 ) 

Q,^(d—k) 


= 0 


and together with (1491) the assertion (l47l) of the theorem follows. 

Because of this result we try further standardization by a d ~ k ~z , it is 


a 2 (d-k)~ 


ipivar G d . k (^) ~ + J 2 - (E G d _ k ( 


First we shall prove that 


lim a 

a—too 



(EG d - k U d) )) 2 

(j 2 (d—k ) 


= 0 , 


(52) 


which follows from Lemma [7] Here the two sums in the Lemma [7] represent 
the two terms in the numerator of (1521) . which in our case can be expressed in 
such a way that each summand is a correlation function of some orientation 
configuration multiplied by I k and by the probability of the occurence of such 
configuration. The convergence rate properties follow from the calculations of 
the denominator of correlation function in Lemma [5j where all parts converge 
at least at exponential rate to some positive value (we omit all zero values). 
The same holds for numerator of correlation function which is calculated in the 
same manner in Lemma [ 6 ] 

In the term J 2 of (1511) we need to consider number of common orientations 
between both functions. Let us have l common orientations e 2 ,..., ez+i not 
including orientation ei of x\, probability of the configuration with this number 
of common orientations is 

d ( d- 1 

d 2 (d-k)-i \i td _ k _i_i^ d _ k _i_i i 2 k + l- d ,+ l 
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and the corresponding correlation function p 2 (d-k)- 1 tends to 2k+l d d+1 in (IMl) . 
Therefore 2 (<i A- t tends to 

, (d — l)(d— k )\ 2 d y^ 1 f d- 2 

k d 2 ^ d ~ k )- 1 \l,d-k-l-l,d-k-l-l, 2 k + l-d 

(d-l)(d-k )\ 2 f d — 2 \ 2 
k d2(d- fc )-i y d _ k _ij > 

where the last equation holds because both sides represent number of possi¬ 
bilities how to choose two sets of distinct d — k — 1 elements from total d — 2 
elements. This corresponds to lim^oo 2(d * fc) . —var G d -k{p i d) ) in BE]). □ 

Remark 5. 77ie standardization of the variance in f f 8 \) corresponds to the 
Poisson case in Subsection 5.1. It remains to investigate whether the central 
limit theorem also holds in our case. 
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